Decoherence and dephasing in strongly driven colliding Bose-Einstein condensates 
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We report on a series of measurements of decoherence and wavepacket dephasing between two 
colliding, strongly coupled, identical Bose-Einstein condensates. We measure, in the strong excita- 
tion regime, a suppression of the mean-field shift, compared to the shift which is observed for a weak 
excitation. This suppression is explained by applying the Gross-Pitaevskii energy functional. By 
selectively counting only the non-decohered fraction in a time of flight image we observe oscillations 
for which both inhomogeneous and Doppler broadening are suppressed, in quantitative agreement 
with a full Gross-Pitaevskii equation simulation. If no post selection is used, the decoherence rate 
due to collisions can be extracted, and is in agreement with the local density average calculated 
rate. 



The dephasing of a momentum excitation in a non- 
uniform Bose-Einstein condensate (BEC) is governed 
by several factors Q, @|. Notably, the inhomoge- 
neous broadening of the Bogoliubov energy shift and the 
Doppler broadening. We denote these processes as de- 
phasing since they involve a reversible, although non- 
linear, evolution of the macroscopic wavefunction. In 
contrast, the decoherence of such an excitation due to 
collisions with the BEC, involves coupling to a quasi- 
continuum of initially unoccupied momentum states 0] . 
This process is driven by the fluctuating ground-state oc- 
cupation of these modes, and is therefore inherently irre- 
versible on timescales larger than the so called " memory 
time" of the system 

In previous works the inhomogeneous lineshape was 
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FIG. 1: Ramsey dephasing calculation. We simulate, by solv- 
ing the cylindrically symmetric GPE [Tlj . a Ramsey type ex- 
periment (solid boxes). We apply a detuned (~ 3 kHz) and 
weak 40 usee Bragg pulse, wait a variable Ramsey delay time, 
and then apply a second pulse which is in-phase with the first. 
A simplified local density averaging is also shown (solid line, 
described in the text) , and is seen to capture the essence of the 
more complete simulation. Both curves decay at roughly the 
expected p| inhomogeneous dephasing rate (3.0 kHz). (lfl | 



seen to agree, for short Bragg coupling times, with both 
the experimental data [g, and with simulations of 
the Gross-Pitaevskii equation (GPE) 0. At longer 
timescales the inhomogeneous lineshape is resolved into 
radial modes |6j. Weak probe Bragg spectroscopy and 
interferometry were also used to characterize the coher- 
ence and spatial correlation function of condensates Q 
and quasi-condensates [9j- In this letter we measure the 
dephasing and decoherence of strongly driven oscillations 
between two identical colliding BECs. By post select- 
ing the non-decohered fraction we observe a strong sup- 
pression of both inhomogeneous and Doppler broadening 
mechanisms, which is in quantitative agreement with the 
results of a GPE simulation. The observed collisional 
decoherence rate is seen to agree with the expected local 
density approximation (LDA) average of the free particle 
collision rate. 

The RMS width of the LDA inhomogeneous lineshape 
for high momentum excitations is given by y8/147/x, 
where fj, is the chemical potential of the BEC The 
RMS Doppler broadening is given by, ^/8/3hk/mR z 0], 
where k is the wavenumber of the excitation, m is the 
mass of the BEC atom and R z is the Thomas-Fermi ra- 
dius of the condensate. For our experimental parameters 
the Doppler dephasing rate is predicted to be 0.6 kHz, 
which is clearly dominated by the expected 3.0 kHz inho- 
mogeneous dephasing rate [ljj- A GPE simulation 
of the final excited population in a weak Bragg/Ramsey 
type experiment (an initial 40 /^sec weak coupling Bragg 
pulse, followed by a delay, and then a second pulse in 
phase with the first), with simulation parameters as de- 
scribed for our experimental system below, is shown in 
Fig. 1 (solid boxes). A simplified LDA of the Ramsey 
signal (solid line), taking only the density inhomogene- 
ity into account, is seen to be in agreement with the 
full simulation. Here the weakly excited state, in the 
Bloch vector picture, precesses around the effective in- 
homogeneous local detuning, and this leads to decaying 
oscillations in the total final excited population. The cal- 
culated inhomogeneous dephasing rate of 3.0 kHz for our 
system agrees roughly with the Ramsey signal both in 
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the simulation and in LDA theory. We conclude that 
weak excitations, at short times, are well understood by 
LDA theory @. 

We now turn to the case of a strong excitation, i.e. 
where population of the zero-momentum condensate is 
completely (or nearly) depleted and transferred coher- 
ently to a travelling condensate in a Rabi-like oscilla- 
tion. For sufficiently large k, the momentum of the 
excitation, we approximate such states by tp(r, z) — 
a4>o(r,z) + be lkz cf>Q(r, z) JL2], where a and b are the am- 
plitudes of the effective two level system and \a\ 2 + 
\b\ 2 = N, with N being the total number of atoms. 
4>o(r, z) is the ground-state wavefunction of the system 
with no excitations. The energy of such a state ip 
can be evaluated by the Gross-Pitaevskii energy func- 
tional 1 3] E = J dV{ £ | (r, z) | 2 + V (r, z) ^ (r, z) \ 2 + 
%\ip(r, z)\ 4 }, where Vo is the external potential and g is 
the mean-field coupling constant. By differentiating the 
energy by the excited state population \b\ 2 , and consid- 
ering the result as a function of |6| 2 , we find a population 
dependant excitation energy. For our experimental sys- 
tem we find a nearly linear decrease in excitation energy 
from the Bogoliubov prediction (weak excitation) at low 
|6| 2 , via the free particle value at \b\ 2 = N/2 down to 
a symmetrically downshifted energy for \b\ 2 « N, where 
the zero-momentum state is simply a weak excitation of 
the travelling condensate. The intuitive picture is that 
the energy per particle required to transfer the conden- 
sate from \b\ 2 = to \b\ 2 = N is simply the free par- 
ticle energy, since the internal interaction energy is not 
changed by such a transformation. At sufficiently large 
Rabi frequencies, that achieve complete population inver- 
sion, the temporary mean field detuning averages to zero 
due to its symmetric nature, and the resonance shifts to 
the free particle value. 

Our experimental apparatus is described in |l4| . 
Briefly, a nearly pure (> 90%) BEC of 1.6(±0.5) x 10 5 
87 Rb atoms in the \F, m/) = 1 2, 2) g round state, is formed 
in a QUIC type magnetic trap [lfij . The trap is cylindri- 
cally symmetric, with radial (r) and axial (z) trapping 
frequencies of 2ir x 226 Hz and 2ir x 26.5 Hz, respectively. 
This corresponds to [i/h = 2.48 kHz. 

We excite the condensate at a well defined wavenum- 
ber using two-photon Bragg transitions The two 

Bragg counter-propagating (along z) beams are locked 
to a Fabri-Perot cavity line, detuned 44 GHz below the 
5S 1 / 2 ,F = 2 — > 5P 3/2 , F' = 3 transition. At this detun- 
ing and at the intensities used here, there are no discern- 
able losses from the condensate due to spontaneous emis- 
sion. The frequency difference du between the two lasers 
is controlled via two acousto-optical modulators. Bragg 
pulses of variable duration and intensity are applied to 
the condensate, controlling the excitation process. 

Following the Bragg pulse, the magnetic trap is rapidly 
turned off, and after a 38 msec of time of flight expansion 
the atomic cloud is imaged by an on-resonance absorp- 
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FIG. 2: Time of flight images of oscillating and colliding 
BECs. (a) A weak perturbation of the BEC (left, in solid 
line ellipse) to the excited population (right, in solid line el- 
lipse). The dashed line marks the region of interest for mea- 
suring the average momentum P to t (b) A n/2 pulse, note the 
strong collisional sphere, (c) Almost a tt pulse, note the weak 
collisional sphere, indicating collisions between the excitation 
and the zero momentum BEC. (d) After further oscillation 
(> 10-7T at 8.6 kHz), the BECs are completely decohered, and 
the Bragg coupling no longer effects the system. 



tion beam, perpendicular to the z-axis. Fig. 2 shows 
the resulting absorption images, for different excitation 
strengths and duration. Fig. 2a shows a perturbative 
excitation with the large cloud at the left corresponding 
to the BEC. A halo of scattered atoms is visible between 
the BEC and the cloud of unscattered outcoupled excita- 
tions to the right. In Fig. 2b we show a tt/2 pulse which 
generates a nearly symmetric excitation. The two con- 
densates collide producing a strong collisional sphere. In 
Fig. 2c we show a nearly complete tt pulse, with a weak 
zero momentum component remaining as an excitation 
of the travelling condensate. We note the weak thermal 
cloud surrounding the origin, which is largely unaffected 
by the Bragg pulse. When we increase the duration time 
of Rabi oscillations (Fig. 2d), the effective two level sys- 
tem is eventually depleted by collisions. 

We measure the response of our system by integrat- 
ing over the elliptical areas shown in Fig. 2a. The solid 
lines are gaussian fit areas to the uncollided atoms. The 
dashed line contains the entire region of interest, includ- 
ing the collisional products. We define the number of 
atoms observed in the excited region inclosed in the el- 
lipse to the right as Nk and those corresponding to the 
initial BEC (left) as Nq. P to t is the total momentum (in 
units of the momentum of a single excitation) along z 
measured inside the dashed ellipse, and normalized by 
the overall number of observed atoms. 
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FIG. 3: Measured Bragg response (solid boxes) spectroscopy 
(N — 0.9 x 10 , and pulse duration of 0.5 msec for these 
measurements). Solid lines are GPE simulations (a) Weak 
excitation spectrum exhibiting the usual mean-field shifted 
Bogoliubov excitation resonance. The dashed line is at the 
theoretical LDA excitation energy (16.15 kHz), (b) Strong 
excitation spectrum, in which a clear suppression of the mean- 
field shift is observed. The dotted line is at the free particle 
excitation energy (15.08 kHz). 



The suppression of the mean-field shift in strongly ex- 
cited condensates is shown in Fig. 3b, where the reso- 
nance is clearly (14.9 ± 0.2 kHz from a gaussian fit) in 
agreement with the free particle value (15.08 kHz, in- 
dicated by the dotted line). This should be compared 
to the resonance (16.05 ± 0.2 kHz from a gaussian fit) 
observed in Fig. 3a for a weak excitation in agreement 
with the expected LDA value of 16.15 kHz (indicated by 
the dashed line). The solid lines are GPE simulations 
of the system, which also confirm the suppression of the 
mean- field shift. 

This suppression of the mean-field shift should cause a 
similar decrease in the inhomogeneous broadening, lead- 
ing to longer coherence times for strongly driven conden- 
sates. We explore this at various driving Rabi frequencies 
(but holding 5uo = 2tt x 15 kHz constant). Fig. 4a and 4b 
show Ptot(t) at the driving frequencies of 1.2 kHz and 8.6 
kHz, respectively. Fig. 4c and 4d show N k /(N Q + N k )(t), 
for the same driving frequencies. 

The dashed lines in Fig. 4a and 4b are exponentially 
decaying oscillations, fitted to the experimental points. 
Here the decay is mainly due to collisions between the 
two condensates. The solid lines in Fig. 4c and 4d are 
GPE simulations with no fit parameters. We observe the 
remarkable result that by post-selecting the non-collided 
fraction we can recover the GPE dephasing behavior, de- 
spite the fact that the GPE totally disregards collisions 
between the excitations and the BEC and considers only 
the evolution of the macroscopic wavefunction. 

The results of this analysis are summarized in Fig. 5. 
Here the solid boxes show the measured dephasing decay 
rate of Nk/(No + Nk)(t). An additional experimental 
point is measured at the Bragg driving frequency of 3.4 
kHz. The solid line represents the theoretical dephasing 
rate obtained by numerically solving the GPE and then 



fitting the time evolution to an oscillating exponential 
decay. 

The empty boxes are the fitted decay rate of the os- 
cillations in P tot (t) . The dashed line is the result of the 
theoretical sum of the GPE dephasing rate calculated 
above (solid line) , plus the free particle LDA collision 
rate (1.66 kHz) |l7| . with no fitting parameters. 
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FIG. 4: Experimental oscillations between two coupled BECs 
(solid boxes) (a-b) average momentum Ptot at 1.2 (a) and 
8.6 (b) kHz oscillations. The oscillations are damped due to 
collisions. The dashed lines in (a) and (b) are exponentially 
decaying oscillation fits to the points, (c-d) Post-selection of 
the uncollided fraction N k /(N + N k ) for the 1.2 (c) and 8.6 
(d) kHz oscillations. The solid lines in (c) and (d) are GPE 
simulations with no fitting parameters. 

There are several points of interest in this figure. 
Firstly, we note that the suppression of inhomogeneous 
effects is again confirmed both in simulation and in ex- 
periment. This can be observed by the minima of the 
simulation curve at 3 kHz driving frequency. The exper- 
imental data agrees with this trend. Where the observed 
decay rates of Ptot and N k /(No + N k ) are 2 and 10 times 
smaller, respectively, than the calculated inhomogeneous 
dephasing rate (as shown in Fig. 1). 

Secondly, the GPE simulation dephasing rate (~ 0.1 
kHz) near 3 kHz driving frequency is significantly slower 
than the Doppler dephasing rate (0.6 kHz). The exper- 
imental points agrees with this trend. This suppression 
can be understood by considering the system in the frame 
of reference of the travelling light potential . In this 
reference frame the condensate is on the edge of the Bril- 
lion zone, and is repeatedly reflected by the potential. 
Due to the strong lattice potential, there is a broad re- 
gion in momentum space for which the group velocity is 
zero. Consequently, the two wavepackets do not separate, 
suppressing the Doppler broadening as well. The reason 
for the increase in the dephasing rate at higher driving 
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frequencies beyond 3 kHz appears to result from the exci- 
tation of other momentum modes along with wavepacket 
spreading in momentum space. 

Thirdly, we note that we can use the measured de- 
cay rates of Ptot{t) and Nk/(N + Nk){t) to isolate and 
estimate the decay due to collisions. We thus measure 
the collisional cross section to be 7.1(±1.8) x 10~ 16 to 2 in 
agreement with the known value of 8.37 x 10 m 
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FIG. 5: Measured decay rates as a function of oscillation fre- 
quency. The solid boxes are the exponential fit decay rates of 
the oscillations in the post-selected Nk/(No + Nk). The solid 
line is the rate of decay fitted to the results of GPE simula- 
tions. The open boxes are the fitted decay rates of average 
momentum Ptot- The dashed line is the sum of the expected 
LDA collision rate (1.66 kHz) with the decay rates from GPE 
simulations (solid line). 

In conclusion, we measure a suppression of the mean- 
field shift for a strong excitation of the BEC. We also 
observe an order of magnitude suppression of the inhomo- 
geneous broadening and Doppler broadening mechanisms 
between strongly driven colliding Bose-Einstein conden- 
sates. Furthermore, we measure a collisional decoherence 
rate in agreement with that expected from previous mea- 
surements. 

In the future we hope to observe a splitting in a Bragg 
probe spectrum from such a strongly driven system. We 
also hope to measure a shift in the energy of the excita- 
tions due to off-resonance Bragg pulses, in analogy with 
the ac Stark shift. This shift may even modify the colli- 
sion rate by shifting the resonance energy sufficiently to 
influence the interaction with the finite width collisional 
quasi-continuum |2c|. 
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